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Abstract 



Using methods of conformal field theory, we conjecture an exact 
form for the probability that n distinct clusters span a large rectangle 
or open cylinder of aspect ratio k, in the limit when k is large. 



The study of the structure of large clusters at the percolation threshold 
continues to pose interesting problems whose solution sheds light on the 
nature of the critical state in general. Recently some attention has been 
paid to incipient spanning clusters (ISCs). These are clusters which connect 
two disjoint segments of the boundary of a macroscopically large region. 
Langlands et al. |l| conjectured that the probability that at least one such 
cluster exists (that is, that the segments are connected) is invariant under 
conformal transformations. This statement was placed in the context of 
conformal field theory in Ref. ||, where an explicit formula was given for 
this crossing probability. (For a review of the status of conformal invariance 
in critical percolation, see Ref. ||.) 

More recently, Aizenman has considered, among other things, the 
probability that there exist n distinct ISCs connecting the two segments.^] 
In the case of a rectangular region, [0, kL] x [0, L], he has proved that the 
probability P(n, k, L) that the strip is traversed (in the direction of length 
kL) by n independent clusters satisfies the bounds 

Ae~ an2k < P{n, k, L) < e~ a ' n2k , (1) 

where a and a' are (different) constants. Note that, on the basis of scale 
invariance at the critical point, P(n, k, L) is expected to have a finite limit 
as L — > oo. 

In this note we extend the arguments of Ref. to determine the exact 
behaviour of the scaling limit of P for large k, namely that 

lim In P(n, k,L) ~ n(n — |) k, (2) 

as k — ► oo for any n. 

An analogous problem may be posed on a open-ended cylinder of circum- 
ference L and length kL. In this case we find, for n > 2, 

lim In P(n, k,L) ~ — — (n 2 - ±) k, (3) 

In passing, we note that, as observed by Aizenman 0], the form of the 
result in (0-0) is not surprising, even though it contradicts what appears to 

1 One may distinguish the probability of exactly n ISCs from that of at least n. However, 
in the limits considered in this note, these will turn out to be asymptotically the same. 
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have been a former consensus that only one such cluster should exist which 
has only recently been challenged and corrected by numerical evidence || [7]]. 
For a recent review see Stauffer ||. Indeed, if one imagines dividing the 
rectangle into two equal rectangles each of size [0, kL] x [0, L/2], and assumes 
that the dominant event will be that of approximately n/2 clusters spanning 
each half, then, up to prefactors, 

P(n,k,L) ~ P((n/2),2k,L/2) 2 . (4) 

Together with the expected exponential dependence on kL at fixed L, this 
leads to the above form. A similar argument then may be made for the 
cylindrical geometry. 

Our argument for the exact coefficients in is based on the well-known 
mapping of bond percolation to the g — *> 1 limit of the g-state Potts model, 
and the understanding of the critical theory of this model through conformal 
field theory. In the Potts model, spins Sj are placed at the sites i of a lattice, 
each taking one of q possible states. The partition function has the form 
Tr Ylij (1 — P+pSsis ), where the product is over all links of the lattice. This 
may be expanded in powers of p/(l — p) so that each term corresponds to a 
particular realisation of bond percolation, weighted by a factor of q for each 
connected cluster. The limit q — > 1 then weights these as in percolation, 
but, as will be seen, it is also often helpful to consider first the case of more 
general q. 

In the rectangular geometry described above, it is useful to express things 
in terms of the transfer matrix e~ H ^ for a strip of width L. The partition 
function for the Potts model with particular boundary conditions at either 
end of a strip of finite length kL then has the form 

(A\e- kL ^ L) \A), (5) 

where \A) is a boundary state corresponding to the boundary conditions cho- 
sen. The symmetry of the Potts model ensures that the degenerate subspaces 
of eigenstates of H(L) may be chosen to transform according to irreducible 
representations of the permutation group S q of q objects. Conformal field the- 
ory also asserts that, in the scaling limit, the states in the low-lying spectrum 
of H(L) transform according to highest weight representations of a Virasoro 
algebra, and their corresponding eigenvalues have the form 7r(x + integer)/L, 
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where x is the highest weight. Thus ([5]) may also be written 

^ e -™nkJ2(A\N)(N\A) e -* Nk , (6) 

R N 

where the first sum is over highest weight representations R, and the second 
over the states \N) is each representation. (This notation is a little corrupt 
because there are in general many states at level N.) 

The simplest non-trivial irreducible representation of S q has dimension 
q — 1, corresponding to a vector (<pi, . . . , ip q ) with J2l=i ¥>a = 0. An example 
is the Potts order parameter ip a = Sg^a—q^ 1 . Out of this other representations 
may be built by taking direct products. For example symmetric tensors (p^ 
with a 7^ b and J2 a Vab = give a representation of dimension (q — l)(q — 2)/2. 
In general we may construct tensors (fab... with n components, none of whose 
indices are equal. Let us denote by R n the Virasoro representation which 
also carries this representation of S q and which has the smallest weight xr. 
Denote this weight by x n . 

Suppose now that we are interested in those configurations in which at 
least n distinct ISCs connect the two ends of the strip. In that case it is 
possible to colour these clusters with n different colours of the Potts model, 
and therefore the states which propagate along the strip must carry at least 
n different colours. If there are fewer than n ISCs, it is not possible to make 
such an assignment. In the limit of large k, then, the partition sum in (|5|) will 
be dominated by those state (s) transforming according to representations R n i 
with n' > n. As we shall argue, the highest weights x n / are monotonically 
increasing in n' . Thus the states with n' = n dominate the sum. 

What is the value of x n l For n = 1 the answer is known, since it cor- 
responds to the scaling dimension of the Potts order parameter near the 
boundary of a semi-infinite system. It was conjectured in Ref. || that this 
corresponds to the operator (1, 3) in the Kac classification || [T0[] , giving, for 
general q, X\ = (m — l)/(m+ 1), where q = 4 cos 2 (7r/(m + 1)), and X\ = | 
for q—1. This conjecture agrees with the known exact result for q = 2 
and numerical work for q = 3 and q—1. Its correctness is also born out 
by the numerical success of the crossing formula of Ref. . We now further 
conjecture that the representations _R n correspond to (1, In + 1) in the Kac 
classification. This is based on the fusion rules for these representations. Ob- 
serve that the composition law for the S q representations under consideration 
is isomorphic to that for addition of spin n in SU(2), which in turn is the 
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same as the fusion rules || for the Kac representations (1, 2n + 1) in confor- 
mal field theory (in non-minimal models corresponding to generic values of 
q). Thus, for example, insertion of two order parameters ip a and ipb near the 
end of the strip will in general give rise to propagating states correspond- 
ing to the tensor representation tp^ (when a ^ b), the vector representation 
(when a = b) and the identity representation. Since these last two corre- 
spond to (1,3) and (1,1) respectively, we may identify the first with (1,5). 
This argument may be generalised straightforwardly to higher values of n. 
Then, according to the Kac formula |§,the highest weight of the (1, 2n + 1) 
representation is x n = n(mn — l)/(m + 1), or x n = n(2n — l)/3 for q = 1. 
This, combined with (§), gives the first result (@), valid as k — > oo at fixed 
n. 

The large n behaviour of (|2|) may also be derived directly from Coulomb 
gas arguments [|11[|. In this approach, the configurations of the critical cluster 



model are mapped onto those of densely packed loops on the surrounding 
lattice. Each loop carries a factor q 1 ^ 2 , which may be traded for local weights 
by considering each loop as corresponding to two oriented loops with vertex 
weights e ±lx according to whether they turn to the right or left at a given 
site, and setting q 1 ^ 2 = 2cos4x- These loop configurations are then mapped 
onto those of a local height model on the dual lattice, with heights 0(r) e |Z, 
and the rule that the height difference between neighbouring dual sites is 
according to the orientation of the corresponding dual bond. This in turn 
is supposed to renormalise onto a Gaussian model with reduced hamiltonian 
(g/47r) J('V4>) 2 d 2 r, where g = 2(2 — 8x/7r), and 2 < g < 4. Free boundary 
conditions on the original Potts model correspond to Dirichlet conditions 
<p = constant in the height model. 

In the strip geometry, the total charge, that is the number of left-oriented 
minus right-oriented loops, is conserved along the strip. Consider the config- 
urations where this charge is 2n. These correspond to cluster configurations 
where at least n distinct clusters traverse the strip, as illustrated in Fig. III. In 
the height model this means that the difference in the heights between the up- 
per and lower edges is fixed to be 2n(n/2) = nn. Neglecting fluctuations, the 
energy of such a configuration is simply (g / 4tt) (mr / L) 2 ■ kL ■ L = (gun 2 /4)k. 
Inserting the q = 1 value g = | gives the leading term in the result @ for 
In P. This calculation works only for large n because (a) the mapping to the 
Gaussian model is valid only in the bulk and not close to the boundary, and 
only for large n are most of the loops far from the boundary; and (b) because 
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Figure 1: Configuration in which two clusters span the strip. The hulls 
of these correspond to a loop configuration with charge 4. Other possible 
non-spanning clusters are not shown. 



the fluctuations are expected to give an 0(1) contribution. 

A similar argument, in this case yielding the exact result at large k, 
may be applied to the cylindrical geometry, where there are periodic bound- 
ary conditions around the strip. Once again, the loop configurations with 
charge 2n correspond to at least n clusters connecting the ends of the cylin- 
ders. Writing <p = rnrv/L + <f)', where v is the coordinate around the cylin- 
der and <$' satisfies periodic boundary conditions, the energy functional is 
(gim 2 /4)k + (g/Air) j(V(f)') 2 d 2 r, where the first term is identical to that for 
the strip with free boundaries. The integral over the fluctuating part then 
gives a contribution ]12| {itcc/fyk to InP, where cq = 1 is the central charge 



of the free scalar field 0'. Putting these together ans setting g = | then 



gives the result in (||]). Note that this is correct only at q — 1: in general it 
should be normalised by the partition function. The finite-size corrections 
to this have the above form, with c = at q — 1. This comes about be- 
cause the Gaussian result cq = 1 is reduced by the effects of loops which can 
wind around the cylinder |13| . These are forbidden when other loops already 
extend along the cylinder, so that no similar reduction occurs in this case. 
In general, the behaviour in (0) should be of the form 2irx$'k, where x$ 



is a bulk exponent JL4J]. In fact, these exponents are the so-called multi-hull 



scaling dimensions discussed by Duplantier and Saleur |DJ. In the plane, 
these determine the power law decay ~ \r± — r<2\ " of the probability 
that two points r\ and T2 lie in the vicinity of the external boundaries, or 
hulls, of n distinct clusters. These are computed in the loop gas in terms of 
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configurations with 2n oriented lines running from r\ to r2, which is precisely 
what we have argued above determines InP for large k. Our result in (||) 
agrees with that of Ref. jnj for these exponents. 



As indicated, @ does not hold for n = 1. This is because a single cluster 
which connects the ends of the cylinder is also allowed to wrap around it: 
this is clearly not allowed for n > 2. For n — 1 the equivalence to the hull 
exponents no longer holds. Instead, we expect for large k that P(l, k, L) is 
asymptotically equal to the probability that the two ends are connected (by 
any number of clusters) and it should behave as exp(— 2irxk), where x = Jj 
is the usual magnetic scaling dimension of the q — 1 Potts model, which 
gives the probability ~ \ri — T2\~ 2x that points r\ and r 2 in the plane are 
connected. Thus for n = 1 on the cylinder, (£J) is replaced by 

\nP(l,k,L) ~ -(5tt/24)A;. (7) 

However, the result in (Q) with n = 1 does have a physical meaning: it is the 
asymptotic probability that the two ends of the cylinder are connected by a 
cluster which does not also wrap around the cylinder. As expected, this is 
much smaller then the unrestricted probability. 

Finally we discuss whether it is possible to compute P(n, k, L) for non- 
asymptotic values of k and n, in the scaling limit L — > oo. This corresponds 
to a generalisation of the calculation of Ref. [|2| , and first involves identifying 
suitable boundary conditions corresponding to the states \A). It is not diffi- 
cult to see that these states should be suitable linear combinations of states 
corresponding to boundary conditions in which each Potts spin is constrained 
to lie in a subset of n states out of the possible q. Let us denote the boundary 
state in which each spin is constrained to take the values a or b or . . . (where 
all the a,b,. . . are different) by \ab . . .). Then P(n,k,L) oc (A n \e- kL " (L) \A n ) , 
where, for example, 

= \a)-\b), (8) 
\A 2 ) = | aft) + \cd) - \ac) - \bd), (9) 

and so on. It may be seen that that these states do indeed transform accord- 
ing to the advertised representations of S q , and so will couple precisely to 
the representations R n , and this will be the dominant coupling in the limit 
k — > oo. However, in order to determine the dependence of P(n, k, L) for 
finite k, in analogy with the argument of Ref. 0, it is necessary to determine 
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the four-point function of boundary condition changing operators which con- 
nect the above boundary conditions to the free boundary conditions along 
the other edges. Unlike the case of Ref. it does not appear that these 
operators, for n > 1, correspond to simple Kac representations. However, it 
may still be possible to conjecture a suitable differential equation or an in- 
tegral representation for this function, as was done recently by Watts |Uj in 
the case of the probability of a simultaneous left-right and up-down crossing 
of the rectangle. A simpler case to consider might be that of clusters which 
span a cylinder of finite length. It is known that in this case the appropri- 
ate matrix elements may be expressed as a linear combination of Virasoro 
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We note that the simple argument given above in (|]) (and the rigorous 
arguments of Aizenman H) provide a simple physical reason why the scaling 
dimensions of composite operators such as those discussed should increase 
like n 2 in two dimensions. The generalisation of Aizenman's argument to d 
dimensions suggests that the rate of increase of — In P(n, k, L) is like n d '^ d ~ 1 ^ . 
However, for d > 2 this quantity is no longer related to scaling dimensions 
by conformal invariance. 

Our conjecture that the relevant scaling dimensions in (0) are the (1, 2n + 
1) operators in the Kac classification is equivalent to a result of Saleur and 
Bauer |Tj| for the spin-n operators in the Bethe ansatz solution of the equiv- 
alent vertex model. These are the 'boundary multi-hull' operators. 

After this work was completed, we saw the paper of Shchur and Kosyakov 



[1£], which reports Monte Carlo measurements of P(n, 1,L) for n = 2 and 
n = 3 on lattices with L up to 64. Their quoted results agree well with 
our predictions in (0, |J), even though the value of k — 1 is not large. In 
particular the ratios of — In P(n, 1, L) between the cases of open boundaries 
(|2p and periodic boundary conditions (|3|) is predicted to be -= — 0.8 for n = 2 
and | ~ 0.857 for n — 3. The corresponding values quoted in Ref. [19 



are 0.808(10) and 0.851(20). This close agreement with the asymptotic form 
may be explained by the observation that the higher eigenstates of H in 
(|5[) give corrections of order e~ 2nk . For n = 1, using (|7]) we find a ratio 
| = 1.6, to be compared with the value 1.5348 extracted from the result 
P(l, 1, L) ~ 0.63665(8) of Hovi and Aharony pOf for the cylinder, and the 
exact result of 0.5 for the square. 

The author thanks M. Aizenman, H. Saleur and L. N. Shchur for useful 
correspondence and discussions. This work was completed while the author 
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